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TECHNICAL NOTE 3456 

PROPAGATION OF A FREE FLAME IN A TURBULENT GAS STREAM 
By William R. Mickelsen and Norman E. Ernstein 


SUMMARY 

Effective flame speeds of free turbulent flames were measured by 
photographic^ ionization-gap^ and photomultiplier tube methods and were 
found to have a statistical distribution attributed to the nature of the 
turbulent field. The effective turbulent flame speeds for the free flame 
^ were less than those previously measured for flames stabilized on nozzle 

burners^ Bunsen burners^ and bluff bodies. The statistical spread of 
the effective turbulent flame speeds was markedly wider in the lean and 
rich fuel -air -ratio regions^ which might be attributed to the greater 
sensitivity of l ami nar flame speed to flame temperature in those regions. 

^ Values calculated from the tiirbulent free-f lame- speed analysis proposed 

by Tucker and the Scurlock-Grover analysis of turbulent flame propaga- 
tion apparently form upper limits for the statistical spread of free- 
f lame-speed data. Hot-wire anemometer measurements of the longitudinal 
velocity fluctuation intensity and longitudinal correlation coefficient 
were made and were employed in the comparison of data and in the theo- 
retical calculation of turbulent flame speed. 


INTRODUCTION 

The high volimietric heat-release rate required in present-day and 
future jet-engine combustors has placed an ever-growing emphasis on t\rr- 
bulent combustion research. The theory of turbulent flame propsigation 
is as. yet in a formative stage^ partly because of the lack of reliable 
experimental methods for determining the effect of turbulence on the 
rate of flame propagation. The purpose of the present investigation was 
to explore a new experimental method believed to approximate more nearly 
a theoretical^ or ideal^ flame model than previous methods. 

The major portion of experimental data on turbulent flame speeds 
has been taken in open flames stabilized on Bunsen burners (refs. 1 to 
3) or on Mache type burners (refs. 4 and 5). A substantial amount of 
experimental data has also been obtained from flames confined in a duct 
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and stabilized on bluff bodies (refs. 6 to 8) . Another experimental 
method consists of measuring the speed of a flame advancing into a tur- 
bulent fuel-air mixture flowing through a tube (ref. 9). In general^ 
the definition of turbulent flame speed used in these methods is that 
given in reference 8: 


Sm = U 

A 


( 1 ) 


(All symbols are defined in appendix A.) Equation (l) has been used 
for local measirrements by the flame -front angle method and for over- all 
flame -speed measurements by the flame-area method. These two methods 
are shown in the following sketches: 



Flame -front angle method Flame -area method 


If 

4 


Turbulent-flame- speed values obtained by these methods may be sig- 
nificantly affected by conditions external to the flame-front proper^ 
such as: 

(1) Large fluctuations in instantaneous flame-front position which 
introduce uncertainty in the determination of mean position in 
long- time-exposure photographs of stabilized flames (ref. 10) 

(2) Existence of "piloting zones" at the rims of Biansen burners_, 
immediately downstream of blviff bodies, or in the boundary layer 
of flame tubes which may affect the burning rate of the flames 
near such zones (refs. 7 and lO) 

(3) Large velocity differences between unburned and biirned gas 
flows which may introduce considerable turbulence in stabilized 
flames (refs. 2, 6, 7, and ll) 

(4) Curvatirre of unb\rmed-gas streamlines at the flame front which 
may introduce substantial error in flame speeds measured in 
stabilized flames by the angle method (refs. 7 and 12) 
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( 5 ) Variations in local flame speeds along the flame envelope which 
may introduce substantial error in average flame speeds measured 
by the flame- area method 

The objective of this investigation was to determine the effect of 
turbulence on the growth of a free flame. It was thought that the free 
flame would not be subject to the complicating conditions present in 
stabilized flames such as those listed. In addition, the free-flame 
growth coiald take place entirely in a homogeneous, isotropic, turbulent 
gj field away from boundary layers and bluff -body wakes, thus providing a 

^ more nearly fundamental flame model. 

The method used in the present investigation is an extension of the 
"soap-bubble" technique (ref. 13) and consists of the observation of the 
growth of a free-flame globule as it is carried downstream in a flowing, 
turbulent, homogeneous mixture of gaseous fuel and air. The flame glob- 
ule was initiated by a single spark, and its growth was recorded by 
three separate methods: high-speed motion picture photography, 

V ionization-gap probes, and photomultiplier tube signals. 

The flame-globule growth was characterized by an effective turbu- 
lent flame speed Sp, which is defined by a mass-balance eqioation similar 

to that derived in reference 13: 


Pa 


( 2 ) 


Equation (2) can “be clarified by the following sketch: 



The rate of change of the globule radius was measured over times less 
than the characteristic time of the turbulence* This means that, on the 
average, the flame front traveled through less than one turbulent eddy 
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during the time of observation. In the past, average turbulent flame 
speeds have generally been measured from photographs having long expo- 
sure times compared with the characteristic time of the turbulence. 

The experimental work was carried out in two phases. The first 
phase was initiated and carried out at the NACA Lewis laboratory and 
is reported in, reference 14. (Because of its limited circulation, con- 
siderable portions of the work reported in ref. 14 are included in this 
report.) Turbulent flame-speed measurements were obtained from high- 
speed motion picture records of the growth of the free-flame globule as 
it was swept downstream in an enclosed tunnel. The turbulent field was 
Introduced by means of a wire grid placed at the tunnel inlet some dis- 
tance upstream of the spark- electrode position. The t\irbulent field 
was characterized by hot-wire anemometer measurements . 

The effective turbulent flame speeds detemined by the photographic 
method varied for runs made under identical fuel-air and mean stream 
conditions. The variation was concluded to be due to the statistical 
nature of the turbulent field. Statistical analysis of groups of 30 or 
more runs made at identical fuel-air and mean stream conditions showed 
that the data for each group approximated a nornnal probability distribu- 
tion. When plotted on the basis of cumulative probability of occiorrence, » 

the effective turbulent flame-speed data showed a consistent increase 
v/ith increasing turbulence intensity. 

The second phase of the investigation was carried out in a free jet 
with a technique whereby the cumulative probability of occurrence of 
turbulent flame speed could be determined from flame groups consisting 
of thousands of separate flame globules. This technique utilized 
ionization-gap probes and counters, which indicated directly the per- 
centage of flames reaching or exceeding any given diameter at a series 
of stations downstream of the spark electrodes. By this method, flame 
speeds were measured over a range of fuel-air weight ratios from 0.053 
to 0.090 and a range of mean stream velocities from 35 to 142 feet per 
second at a stream static pressure of 1 atmosphere and a stream static 
temperature of 85 P. The stream turbulence was varied by the use of 
three interchangeable grids, and the parameters of the turbulence were 
measured with constant- temperature hot-wire-anemometer instr\mientatlon. 

The Investigation was carried out at the Lewis Flight Propulsion 
Laboratory as a part of the combustion research program. 
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APPARATUS 

Air-Ducting^ Wind-Tunnel^ and Air-Flov Instrumentation 

The air-ducting^ wind-timnel^ and air-flow Instrumentation used in 
the first phase of the investigation and reported in reference 14 is 
shown in figure l(a) . Alterations made to the apparatus for the second 
phase of the investigation are shown in figure l(h). The alterations 
essentially converted the system from an enclosed tunnel to a free jet. 

Air supplied from the laboratory air facility was metered with a 
standard variable-area orifice and controlled with the butterfly and 
gate valves shown in figure 1. Static pressures and temperatures were 
measured with conventional manometers and thermocouples. 

The ducting approach to the tunnel was designed from principles 
used in conventional low- turbulence wind tunnels. The transition sec- 
tion was followed by a honeycomb and turbulence -damping screen in the 
calming section to reduce the turbulence present in the air-supply flow. 
The entrance to the enclosed tunnel (or free jet) was a nozzle with a 
contraction ratio of 20^ and provisions were made at the tunnel throat 
» for the insertion of turbulence -producing grids o Three grids were used 

for the investigation^ all having a mesh-to-wire diameter ratio of 5. 

The grid wire diameters were 0.0313^ 0.063^ and 0.125 inch. 

The enclosed- tunnel installation was equipped with an adjustable 
plug and nozzle by which critical flow could be maintained at the tunnel 
exit^ so that the tunnel flow field could be isolated from pressure dis- 
turbances in the laboratory exhaust facility. The free-jet installation 
was equipped with an exhaust hood into which diluting room air could be 
drawn in addition to the free-jet flow. 


Fuel System 

The propane fuel used throughout the investigation had the follow- 
ing composition by liquid volume: 


Propane^ percent 97.76 

Ethane^ percent 1.72 

Isobutane ^ percent .51 


The propane fuel system is shown in figure 2 and was identical for both 
the enclosed- tunnel and free-jet installations. Propane was supplied 
from a laboratory facility and metered with standard rotameters. Pres- 
sure ^ temperature, and flow controls are shown in figure 2. The propane 
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was injected into the air stream through a grid having 60 equally spaced 
0.040- inch- diameter holes. The injection grid was installed at a station 
upstream of the calming chamber as shown in figure 1. 


Ignition System 

The flame-globule ignition system consisted of a pair of spark 
electrodes protniding into the tunnel stream from opposing walls and an 
electrical energy soiirce. The spairk electrodes used in the enclosed 
t\mnel were carefvilly streamlined to a fineness ratio of 5 with a thick- 
ness tapering from 0.019 inch at the tunnel center to 0.038 inch at the 
tunnel wall. The spark electrodes used in the free- jet installation 
were of circular cross section with a diameter tapering from 0.015 inch 
at the stub tunnel center to 0.075 inch at the stub t unn el wall. The 
spark-gap spacing was 0.015 and 0.030 inch for the enclosed- tunnel and 
free-jet installations, respectively. The spark- electrode positions 
for the two installations are shown in figure 1. 

The spark-energy so\rrce was of the capacitive type and produced 
single sparks synchronized with the camera operation in the enclosed- 
tunnel installation or successive sparks at three per second in the 
free-jet installation. The spaxk energy could be varied by changing 
capacitors . 


Flow-Field Instrumentation 

In general, the flow-field instrumentation was adaptable to point 
measurements in either the enclosed- tunnel or free-jet streams. Mean- 
stream-velocity measiirements were made with a conventional pitot-static 
probe and mlcrcmanometer. Fuel-air ratio was meas\ared directly by a 
mixture analyzer. Stream samples were drawn through a sharp-edged 
0.125-inch- inside-diameter sampling probe to the analyzer, which indi- 
cated fuel-air ratio directly by the thermal- conductivity-bridge method. 

Turbulence measurements were made with constant -temperature hot- 
wire -anemometer equipment, which is described in appendix B. Sound- 
pressure levels were measured with a conventional microphone and a 
sound-pressure-level meter, also described in appendix B. 


Free-Flame-Growth Measurement 

Photographic instrumentation . - The growth of the free -flame glob\ile 
was recorded photograpMcally in the enclosed- tunnel installation with 
the optical system shown in figure 3. The photographic records were 
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taken with a 16 -millimeter Fastax camera operated at approximately 3600 
frames per second and located at the end point of the parallel-heam 
schlieren system. The enclosed- tunnel windows were strlation-free plate 
glass. The light source was a mercury- vapor BH-6 lamp operated with a 
high-voltage direct -current power supply. 

Ionization-gap ins t lamentation . - The gro\>rth of the free-flame 
globule was measured in the free-jet installation with an ionization- 
gap probe^ a photomultiplier flame sensing unit^ and two flame counters. 
The physical orientation of these instnamentation components is shown 
in figures 4(a) and (b). The photomultiplier unit and its accompanying 
flame counter recorded the total number of flames that occurred during 
any particular run. The ionization -gap probe and its accompanying flame 
counter recorded the number of flames that grew to a radius of r or 
greater at the axial station x during any particular run. 

The photomultiplier unit consisted of a 931A photomultiplier vacuum 
tube with a fixed-plate voltage supply. The ionization-gap probe is 
shown in detail in figures 4(c) and (d) . The brass fairing along the 
probe tip prevented flame seating on the probe both by flame quenching 
and by elimination of the recirculation zone. Because of the convoluted 
form of each flame globule^ the voltage signal from the ionization gap 
usually contained multiple peaks ^ which caused spurious counts in the 
flame counter. In order to smooth these peaks ^ resistors^ capacitors, 
and diodes were included in the 360-volt direct- current battery circuit 
for the ionization gap. 

In order to prevent spurious counts arising from small fluctuations 
in the voltage signal from both the photomultiplier unit and the 
ionization -gap probe, a long time constant was included in the flame- 
counter design. These flame counters employed thyratr on- controlled re- 
lays and solenoids which actuated a mechanical counter. They were de- 
signed by C. C. Conger of the Lewis laboratory^ 

Photomultiplier instrumentation . - Simultaneous measurements of the 
spark light intensity and ensuing free -flame -globule growth were made 
in the free-jet installation with the photomultiplier instrumentation 
shown in figure 5. The photomultiplier unit consisted of a 93IA photo- 
multiplier tube with an adjustable -plate voltage supply. Light emitted 
from the spark and ensuing flame passed through the narrow slit to the 
photomultiplier imit. The voltage signal from the photomultiplier unit 
was recorded on 35-millimeter film from an oscilloscope trace. The 
photomultiplier instrimientation and the assumption required to relate 
output voltage to flame speed are discussed in appendix C. 



8 


NACA TW 3456 


PROCEDURE 

The test schedule was as follows: 


Grid size. 

Stream 

Propane -air 

Distance 

in. 


velocity^ 

weight ratio 

downstream 





of grid^ 

Wire 

Mesh 

ft/sec 



diam. 




in. 

Free jet 

0.125 

0.625 

35 to 105 

0.07 

15.2 to 16.2 

.125 

.625 

70 

.056 to .090 

15.2 to 16.2 

.063 

.313 

70 

.058 to .084 

15.2 to 16.2 

.0313 

.156 

70 

.053 to .084 

15.2 to 16.2 

Enclosed tunnel 

0.125 

0.625 

35 to 140 

0.045 (approx.) 

6.3 to 15.3 


Enclosed- Tunnel Installation 

Room air was supplied to the enclosed tunnel and exhausted at 
critical flow condition through an exit choke by means of the altitude 
exhaust facility. Timnel static pressxires and temperatures were meas- 
ured by manometers and thermocouple rakes as shown in figure l(a) . Fuel 
was injected into the air stream and metered through a standard rotam- 
eter^ after which dilution air was admitted downstream of the choke. 

The tunnel stream velocity was then measured with a pitot-static probe 
and micromanometer. 

The sequential schlieren photographs of the propagation of the free- 
flame globule in the fuel-air mixture were obtained by actuating a 
single switch which synchronized the following operations : fuel shut-off 

(producing a semi-infinite fuel slug); camera start; ignition by a 
single spark; and camera stop after an appropriate interval. 


Free-Jet Installation 

An exhaust air-flow rate of approximately 500 pounds per minute 
was initiated, and dilution air was admitted to the hood. Laboratory 
pressurized air and propane were admitted into the stub tunnel, after 
which the fuel-air ratio was adjusted to the desired value by meains of 
a probe and the NACA mixture analyzer; the stream velocity was deter- 
mined from pitot-tube measurements. 
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The ionization -gap probe was positioned at a station x downstream 
of the spark electrodes and at a radial distance r from the tunnel 
centerline. A series of flame globules was initiated by the ignition 
system, firing at the rate of three sparks per second. The number of 
flame globules Np intercepting the ionization-gap probe and the total 

number of flame globules ignited were read from instrimientation 

shown in figure 4(b). 


FLOW-FIELD MEASUEEMENTS 

In order to check the uniformity of the flow field, velocity pro- 
files were measured in both the enclosed- tunnel and free- jet installa- 
tions. The mean stream velocity for both installations was essentially 
constant over the core of interest, as shown by figure 6. Similar 
measurements were made to check the uniformity of the propane -air ratio 
in both installations. The propane-air ratio was found to be constant 
over the major portion of the enclosed -tunnel flow field, and typical 
profiles measured in the free -jet installation are shown in figure 7. 
cu 

In order to establish the characteristics of the tiirbulent field, 

^ extensive measurements of the turbulence intensity the spectrum 

function F, and the correlation coefficient f were made in both the 
enclosed- tunnel and the free- jet installations. The definitions of 
these quantities and the instrumentation used to measure them are dis- 
cussed in appendix B. 


Enclosed-Tunnel Turbulence Measurements 

Longitudinal turbulence intensity measurements made in th^ 
enclosed- tunnel installation are shown in figure 8. The data shown in 
the figure agree well with the von Karman data quoted in reference 15 
and show a substantial decrease in turbulence intensity from the elec- 
trode position (x/d = 50) to the end of the flame-growth observation 
field (x/d = 122). Longitudinal intensity profiles measured in the 
tunnel center plane perpendicular to the observation plane were essen- 
tially flat over the core of interest. Turbulent energy spectra were 
measured in the enclosed- tunnel installation as described in appendix B. 
The scale of turbulence was estimated by comparison of experimental data 
with a grid of spectrum density curves based on an e:monential form of 
the longitudinal correlation coefficient f = exp (^/l^ ) . A typical 

measured spectnmi is shown in figure 9, and the estimated turbulence 
scales are listed in table I. 
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Free-Jet Turbulence Measurements 

Conversion of tbe enclosed tunnel to a free jet resulted in an in- 
tense sound field originating in the dilution air and exhaust hood and 
piping. As described in reference 16, velocity fluctmtions associated 
with sound waves contribute to the hot-wire anemometer measurements 
through the relation 



with the assumption that no correlation exists between the turbvilence 
and sound velocity fluctuations. In order to determine the applica- 
bility of equation (3) in the free jet, longitudinal velocity fluctua- 
tion intensities were measured both with and without the intense sound 
field. These measurements are shown in figure 10 for each of the three 
turbulence grids used in the investigation. The velocity fluctuation 
intensity was definitely increased with the presence of the sound field 
for each of the three turbulence-generating grids. 

The sound-pressure level in the free jet was measured with the 
instrumentation described in appendix B and was found to be 136 decibels, 

which corresponds to a soimd- velocity fluctuation of 0.96 foot per 

second. Values of the sound-velocity fluctuation /\^u| were also cal- 
culated from the data of figure 10 by the use of equation (3) . These 
values of are compared in table II with the value of cor- 

responding to the sound— pressure level of 136 decibels j good agreement 
is shown in every case. 

The kinetic-energy spectrum of the sound field alone is shown in 
figure llj it has a generally continuous form with peaks, or periodici- 
ties, at 57, 110, 140, 170, 230 to 285, 600 to 870, 1050, 1500, 2400, 
3000, 4000 to 5600, and 7000 to 9000 cycles per second. Periodicities 
observed in the sound spectrum were generally found again in the longi- 
tudinal velocity fluctuation energy spectra measured in the free- jet 
installation, as shown by the typical spectnmi in figure 12. 

Longitudinal double-velocity correlation coefficients f measured 
in the free-jet installation are shown in figure 13 for each of the 
three turb\ilence-producing grids. The irregular form of the correlation- 
coefficient curves is attributed to the periodic nature of a considerable 
portion of the velocity fluctuations. The correlation coefficients were 
measured by the special methods described in appendix B and are defined 
by the relation 
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( 4 ) 


where the symbols are defined by the following vector diagram: 




U 


OQ 

cn 

CM 


The longitudinal turbulence scale is usually defined by the 
relation 


Because of the periodic nature of the combined turbulence and sound 
fields^ the longitudinal correlation coefficient did not approach 
zero within the range of measurement. Therefore^ equation (5) could 
not be used to define the turbulence scale. Since the growth of the 
flame globule was observed over times comparable with C = 0.1 inch^ 
only the first portions of the correlation curves shown in figure 13 
were considered relevant to the experiment. As a part of the calcula- 
tions involved in the Scurlock-Grover analysis^ values of the lateral 
scale of turbulence Lg were calculated as described in appendix D. 
These calculated values of Lg are listed in the table included in 
the RESULTS section of the text. 


The relation between flame speed and free-flame globule- expansion 
rate is derived in reference 13 for a constant-pressure laminar free 
flame propagating in a quiescent fuel-air mixture. This relation was 
obtained by equating the mass flow of the unbumed gases entering a 
flame front of infinitesmal thickness to the mass flow of the burned 
gases leaving the flame front. In a physical sense ^ the following equa- 
tion accounts for both the flame-front motion due to flame propagation 
and the motion due to thermal expansion of the burned gases: 



( 5 ) 


DETERMINATION OF FLAME SPEED 


Fundamental Considerations 



Pf dr 


( 2 ) 


12 


NACA TN 3456 


The following modification of this equation to account for a finite 
flame-front thickness 5 has been presented in reference 14: 


For this investigation^ the thickness of the flame front was con- 
sidered negligible compared with the radius r so that eq,\mtion (2) 
was used to analyze the data. The values of Pg., computed from Rg 
and Tg, were obtained from measurements of tunnel approach-stream 

flow. The values of computed from Rp and Tj, for each fuel-air 

ratio, were taken from references 17 and 18, respectively. The value 
of Pp/Pg varied from 0.12 to 0.14 in the present investigation. 

A mmiber of assumptions are necessary to the definition of effec- 
tive turbiilent flame speed as given by equation (2) . The definition 
req.\iires that 

(1) The free-flame globule grow uniformly in all radial directions 

(2) The combustion proceed in a constant pressure field 

(3) The flame globule be swept downstream at the stream velocity 

(4) No circimiferentlal flow exist about the flame globule 

(5) Variations in spark energy not affect the globule propagation 


The validity of the assumptions was established by special tests 
described later in this section. 


Photographic data obtained in the enclosed tunnel were analyzed by 
projecting the 16-milllmeter motion plctxore film record of each flame- 
globule expansion onto a film viewer. A typical photographic record is 
shown in figure 14(a) . The outline of each globvile was traced in order 
to permit planimetering its area and calculating a mean radius. The 
mean radius was plotted against time for each flame globiile of the se- 
quence, as shown in figure 14(b). The slope of a straight line faired 
through such a series of experimental points is dr/dt. The effective 
turbulent flame speed was then calculated by multiplying dr/dt by 
the expansion factor Pp/Pg^ shown in equation (2). 



( 6 ) 


rate 


Photographic Method 
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From sets of 30 or more flame speeds obtained in this manner for 
identical fuel-air-flow conditions, the cumulative probability of oc- 
ciirrence of flames having these flame speeds was calculated. In order 
to Illustrate the statistical nature of these data, the turbulent flame 
speeds of a typical set of data have been plotted against cumulative 
probability of occurrence in figure 14(c). Each data point represents 
the effective flame speed of a single free flame. Cumulative proba- 
bility of occvirrence is defined as the percentage of flames having flame 
speeds greater than a given value. 


Ionization-Gap Method 

In order to obtain a set of turbulent flame speeds for a particular 
condition of turbulent approach- stream Intensity and fuel-air ratio in 
the free- jet installation, the ionization-gap probe was positioned con- 
secutively at two stations x = 15.2 and 16.2 Inches, as outlined dla- 
grammatlcally in figure 15. For each chosen position r, at least 
1000 flames N-^ were ignited. The ionization-gap counter recorded the 

number of flames Np intercepting the probe. The ratio Np/N^ was the 

cumulative probability of a flame globule having a radius at least as 
great as that indicated by the probe setting r. The ratio Np/N^ was 

then plotted on probability coordinates against r and a curve was 
faired through each of the two sets of data as shown in figxore 16. The 
difference Ar was obtained from the cuorves for the emulative proba- 
bilities of 2, 20, 50, 80, and 98 percent. Equation (2) was then em- 
ployed to determine the flame speeds corresponding to the five proba- 
bility percentages. 


Photomultiplier Method 

The photomultiplier method was used to obtain values proportional 
to the effective flame speed. The method essentially consists of meas- 
uring the free-flame globule -growth rate along the axial diameter of 
the globule. The method is described in the APPARATUS section and in 
appendix C. The instrumentation is shown in figure 5. 


Verification of Experimental Method 

The discussion of equation (2) in a preceding section included a 
nmber of assmptions necessary to the derivation of the equation. 
Since equation (2) essentially expresses the experimental method used 
in the investigation, the assmptions were verified by tests that are 
described in the following paragraphs. 
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Uniformity of globule growth . - Photographic data of the flame- 
globule growth in the enclosed tunnel were analyzed to determine 
whether the globules had preferential growth directions. This analysis 
was made by measuring globule diameters (see fig. 17) at four angles 
of rotation from the tunnel axis. While each globule had different 
growth rates along the four diameters, no consistent relation between 
growth rate and direction was found. This observation supports the use 
of a mean radius in the photographic method of flame-speed measiorement, 
the use of a single radial direction in the ionization-gap method, and 
the use of a single diameter in the photomultiplier method of flame - 
growth observation. 

Pressure pxilsatlons. - In order to determine whether the fluctua- 
tions in flame -globule -growth rates, illustrated by figure 14.(b), were 
due to pressure fluctuations in the enclosed tunnel, data as shown in 
figvire 17 were inspected for possible phase relations between the 
fluctuations in growth rate along the four diametral directions. Since 
no frequency-phase relation was found from these data, it was concluded 
that pressure fluctuations were not responsible for the fluctuations in 
globule -growth rate. Additional measurements with a wall-mounted pres- 
sure pickup indicated the absence of significant pressure pulsations in V 

the enclosed tiinnel. 

Axial velocity of flame globule . - The effect of buoyancy of the 
hot flame globule on its axial velocity was determined by measuring the 
distance from the spark electrodes to the globule center on each frame 
of photographic sequences such as shown in figure 14(a). The axial 
velocities calculated from such measirrements agreed closely in all 
cases with the measxared stream velocity, which indicated that the 
buoyancy effect was negligible. 

Unburned- gas motion . - As shown by equation (2), the rate of 
thermal expansion of the flame globule is from 7 to 8 times as great 
as the effective turbulent flame speed. In order to determine whether 
the unbvirned-gas motion caused by this thermal expansion was pixrely 
radial (from the globiile center) , simultaneous oscillograph records 
were made in the free- jet Installation with the ionization-gap probe 
and a shielded hot-wire anemometer probe as shown by the typical oscil- 
lograph and probe position diagram in figure 18. The ionization- gap 
trace indicated the flame -globiile passage, and the anemometer trace was 
proportional to the Instantaneous stream velocity. Inspection of 60 
such records indicated no axial change in mean stream velocity and there- 
fore no circumferential flow of unburned gases about the sphere. The 
marked decrease in amplitude of velocity fluctuations at the time of 
flame-globule passage is interpreted as being due to the radial transla- 
tion of the free- jet stream core by the thermal expansion of the flame 
globule . 


3637 


NACA TN 3456 


15 


Effect of spark energy . - In order to detennlne the effect of 
spark energy on the flame -globule growth^ two series of r\ms were made 
in the free-jet installation under identical stream and fuel conditions 
but with spark source energies differing by a factor of 2.5. As shown 
in figure 19^ cumulative probabilities of flame -globiile sizes were 
measured at a series of axial positions with the ionization -gap instru- 
mentation for the two spark source capacitors of 0.0004 and 0.0010 
microfarads. The data of figure 19 show an increase in initial flame- 
globiile size, but no change in turbulent flame speed between the axial 
positions x = 15.2 and 16.2 inches. 

Effect of spark light intensity . - Additional information on the 
relation between the spark characteristics and the flame-globule growth 
was obtained with the photomultiplier instrumentation shown in figure 
5. With this instmmentation, oscillographs were obtained in a free- 
jet installation showing spark light Intensity and flame light intensity 
as a function of time for 151 consecutive flame globules at constant 
mean stream and fuel conditions. Typical oscillographs obtained by 
this method are shown in figure 20. As shown in appendix C, the por- 
tion of the trace due to the spark light intensity was proportional to 
spark cxirrent, and statistical analysis of the magnitudes of the peak 
spark light intensity indicated that this parameter approximated a 
normal probability distribution. The flame speed of the globule was 
taken as proportional to the slope of the portion of the oscillograph 
trace due to flame light intensity, as discussed in appendix C. 

The group of 151 consecutive runs are shown in figure 21, in which 
peak spark light intensity is plotted against flame speed for each par- 
ticular flame globule. The statistical correlation coefficient for this 
data was determined with Pearson's product moment formula (ref. 19) and 
was found to have an absolute value of |0.17| for a possible range of 
zero to 1 1.0 I . The low value of the statistical correlation "coeffi- 
cient is interpreted to indicate that variations in spark cxirrent had 
no effect on txirbulent flame speed. 

Additional verification of statistical variation of turbulent 
flame propagation . - The cxmixilative probability of occxrrrence of slopes 
measured from the series of 151 photomxiltiplier oscillographs is shown 
in figure 22. The statistical distribution of flame speeds measxored by 
this method is interpreted as additional proof of the randomness of the 
phenomena . 

Validity of common origin assumption . - The ionization-gap 
method, described earlier in this section, assumes that all flame 
globules have a common origin point, as shown in figure 15. In 
order to determine the error introduced by this assumption, "origin 
points" were measured from the series of 151 photomultiplier traces by 
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extrapolating the slope of the flame light intensity trace to the time 
axis^ as shown in figure 20. The origin points showed a statistical 
deviation in position, and the theoretical error introduced into flame 
speeds determined by the ionization-gap method was calculated to be 
±2.5 percent or less, within one standard deviation of origin position. 
The error in flame speed due to shift in origin point is small compared 
with the deviations in flame speed actually measured by the photographic 
(fig. 14(c)), the ionization-gap, and the photomultiplier tube methods 
(fig. 22). 


Precision of Data 

Ionization- gap measiirements . - Each data point plotted as cumula- 
tive probability against flame radius, such as shown in figure 16, 
represents information from at least 1000 flames. As shown in refer- 
ence 20, samples consisting of 1000 flames would have the following 
accuracy in the cumulative probability value; 


Cumulative 

Possible 

probability. 

error, 

percent 

percent 

50 

±5 

20 

±7 

2 

±29 


The possible errors listed correspond to 1.28 standard deviations of 
the error distribution and therefore should hold 80 percent of the 
time. The general consistency of the data throughout the lower prob- 
ability region such as shown in figure 16 seems to justify greater 
confidence in the low cumulative probability range than that shown in 
the preceding table . 

The reliability of the curves faired through the plotted experi- 
mental data, such as presented in figure 16, may be estimated from the 
deviation of data points from the fitted curves. ■ 

The greatest divergence was found in results for the 0.125-inch- 
dlameter screen. For the 98-percent probability level, the maximum 
data scatter from a faired curve was 20 percent; for the 50-percent 
level, 7.8 percent; suid for the 2 -percent level, 3 percent. The maxi- 
mum scatter for the 0.063-inch-diameter screen was 5.3 percent at 98- 
percent probability, and for the 0.0313- inch- diameter screen, it was 
13 percent at the 98-percent level. 
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As discussed earlier in this section, the value of flame tem- 
perature used in the equation for turbulent flame speed was a cal- 
culated theoretical value. It may be of interest to note that flame 
temperature has a much greater effect on the laminar flame speed than 
on the calculation leading to effective turbulent flame speed. As 
shown by reference 21, a decrease of theoretical adiabatic flame 
temperature from 4000^ to 3900^ R would result in a laminar flame 
speed decrease of 43 percent. In contrast, this change in flame 
temperature would cause only a 2.5-percent decrease in the turbulent 
flame speed calculated from the equation 


_ PF dr 
^ Pa 


( 2 ) 


It appears, then, for the purpose of this investigation, that the ef- 
fect of a small error in actiml flame temperatiore is negligible in 
calciilating the density ratio. 

Other sources of inaccuracy may be found in the instrumentation 
and ecLuipment. The positioning of the ionization-gap probe involved 
an error not exceeding 0.5 percent. Measurement of the fuel-air ratio 
was accurate within 2 percent. Flow measurements were made with a pre- 
cision of 3 percent. For the measurements characterizing the turbu- 
lence parameters , an estimated value of ±5 percent is suggested for the 
probable limit of error. 


Photomultiplier method . - A treatment of the sources of error 
present in flame-speed measurement by the photcmioltiplier method is 
presented in appendix C. 


Photographic method . - In addition to the errors introduced by the 
equipment and instrumentation, sources of error peculiar to the photo- 
graphic method must be considered. Such factors as adjustment of the 
projected image to one-to-one correspondence, timing, optical resolu- 
tion, and scale factor are all considered significant and have been 
thoroughly treated in reference 22, which assigns a relative maximum 
error of 10 percent to the photographic -reduction technique. 


RESULTS 

Effective turbulent flame speeds were measured in the enclosed- 
tunnel and free-jet installations over a range of propane-air and tiirbu- 
lence conditions as shown in the following table: 
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Grid size, 
in. 

Stream 

velocity^ 

ft/sec 

Longitudinal 

velocity 

fluct\iation 

intensity, 

aIv^, 

ft/sec 

Lateral 

turbu- 

lence 

scale, 

in. 

Propane - 
air 

weight 

ratio 

Distance 
downstream 
of grid, 

in. 

Diam. 

Mesh 

Free jet 

0.125 

0.625 

35 to 105 

1.18 to 2.60 

0.0645 

0.07 

15.2 to 16.2 

.125 

.625 

70 

1.69 

.0645 

.056 to 

15.2 to 16.2 






.09 


.063 

.313 

70 

1.34 

.0587 

.058 to 

15.2 to 16.2 






.084 


.0313 

.156 

70 

1.17 

.0437 

.053 to 

15.2 to 16.2 






.084 


Enclosed tunnel 

0.125 

0.625 

35 to 140 

0.81 to 3.51 

0.08 to 

0.045 

6.3 to 15.3 





0.13 

(approx.) 



The lateral turbulence scales L„ for the enclosed tunnel were obtained 
from the values of Lp shown in table I by using eq.uation (B7) of ap- 
pendix B. The values of Lg for the free jet were obtained by fitting 

an exponential curve to the lateral correlation coefficient g as de- 
scribed in appendix D. Both methods for evaluation of Lg are based 

on the assumption of isotropy, which was not proven because of the lack 
of turbulence measurements in the lateral direction. 


The effective turbvilent flame speeds, plotted in figure 23, were 
meas\nred in the free- jet installation at a fixed propane-air ratio of 
0.07 over a range of mean stream velocities with the 0.125-inch wire 
diameter by 0.625-inch mesh turbulence-generating grid. The data are 
plotted as a function of longitudinal velocity fluctuation intensity 

with emulative probability as a parameter. Since turbulence 
measurements were made at only one stream velocity, and since the 
sound- pressure level was independent of stream velocity , the data of 
table II and the following equation were employed to calculate the 
longitudinal velocity fluctuation intensity: 




+ 


(V) 
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The turbulence intensity was assumed constant throughout the 

range of stream velocity U (see ref. 15), and u2 corresponded to the 

s 

measured so\md-pressure level of 136 decibels. 


Effective turbulent flame speeds measured in the enclosed- tunnel 
installation at a fixed propane-air ratio of approximately 0.045 are 
plotted in figure J4 as a function of the longitudinal velocity fluctu- 


ation intensity Vu“^ with cumulative probability as a parameter. The 
data show a consistent increase of effective turbulent flame speed with 
velocity fluctuation intensity . These data were obtained downstream of 
the 0.125-inch wire diameter by 0.625-inch mesh grid over a range of 
stream velocities from 30 to 140 feet per second. As shown by the data 
of figure 8, the decay of the turbulence intensity was quite severe 
over tte r^ge of interest (x/d = 50 to 122). A turbulence intensity 

of == 0.027 was taken as a representative value and assimied 

constant to obtain values of u2 over the stream-velocity range . 


The effects of propane-air ratio and turbulence scale on the propa- 
gation rate of the free flame were investigated in the free-jet instal- 
lation by measiiring the turbulent flame speed downstream of the three 
turbulence-generating grids over a range of propane-air ratio from 
0.053 to 0.090. (The variations in tiorbulence scale are listed in the 
table at the beginning of this section.) These data are shown in fig- 
ure 25 along with a dashed line representing laminar flame-speed data 
reported in reference 21. Txarbulent flame speeds from the enclosed 
tunnel are included in figure 25(a) and were evaluated by interpolat- 
ing the data of figure 24 at a velocity fluctiiation intensity equal to 
that of the free-jet data in figure 25(a). 


DISCUSSION 


Measurements of Present Investigation 


A discussion of the effect of turbulence on flame propagation is 
generally conceded to rest on the principles of tvirbulent motion set 
forth by Taylor in reference 23. These principles are summarized in 
the equation relating the mean displacement X of a fluid particle 
from its original position to two parameters of the turbulent field; 


D 


1 dX^ 

2 dt 



0 


( 8 ) 
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where D is the turbulent diffusion coefficient in the x-direction, 
X is the particle displacement in the x-direction^ and 0t is the 
Lagrangian double -velocity correlation coefficient defined by the 
relation 



It is of particular importance to note that the turbiilence pajrameters 
in equation (8) are both mean values, and that instantaneous, local 
particle displacements depend on the statistical deviations about those 
mean values. It is suggested that the fluid particle displacements have 
a statlsticaJ. distribution about the mean much the same as the statisti- 
cal variation of velocity fluctuations in the turbulence. Such varia- 
tions in particle displacements woxild be particularly noticeable when 
observed over times less than the characteristic time of the t\arbulence, 
as is the case in the present investigation. 

A statistical variation in the particle displacement X would 
probably cause variations in the speed of flames where flame wrinkling 
was controlling; a statistical distribution of the turbulent diffusion 
coefficient D would also probably cause vaxiations in the flame speed 
where the diffusion process affects turbulent flame propagation. As a 
consequence of this statistical distribution of turbulence velocity and 
correlation coefficients, particular free flames might encounter condi- 
tions especially favorable and therefore have high propagation rates, 
while others might be subjected to conditions that would reduce their 
rate of growth. This picture of statistical variation is consistent 
with the data of the present investigation shown in figures 23 to 25 
and also agrees with obse* vations of stabilized flames that show a wide 
brush of flame fronts. 

The effect of flame temperature on laminar flame speed has been 
reported in references 21 and 24, wherein it is shown that small changes 
in flame temperature result in large changes in laminar flame speed. 

In particular, reference 24 shows that this dependence of laminar flame 
speed on flame temperature is of greater magnitude in the lean and rich 
regions. The lower portions of the flame-speed probability band in 
figure 25 generally show decreases in the ratio of the tixrbulent to 
the laminar flame speed Sp/S^ In the lean and rich fuel-air regions. 

This decrease is consistent with the effect of flame temperature if the 
lower portion of the probability band is considered to represent flames 
that have encountered turbiilent conditions sufficiently violent to de- 
crease the flame temperature by diffusive action. 
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Comparison of figures 25(a) ^ (b), and (c) shows that greater por- 
tions of the flame -speed band fall below the laminar flame speed as the 
grid size becomes smaller. As shown by figure 13^ decreases in the grid 
size reduce the mean tiirbulent eddy size. Since the effect of the t\nr- 
bulence is tho\ight to be increasingly diffusive in character as the tur- 
bulence scale decreases, this trend might be due to an increasing propor- 
tion of flames having temperatures less than the adiabatic value as the 
turbulence scale becomes smaller. 


A comparison of the data of the present investigation with those 
of references 4, 8, 9, and 22 is presented in figure 26. In order to 
obtain a basis for comparison, the data are plotted with the ratio of 
turbulent to laminar flame speed St/Sl against longitudinal velocity 

fluctuation intensity The turbulent flame speeds measured in a 

Mache (nozzle) type burner (ref. 4) , a Bunsen burner (ref. 8), and a 
flame tube (ref. 9) fall above those measured by the free-flame method 
in references 14 and 22 and in the present investigation. It should be 
noted that the Sji/Sl data from the present investigation are plotted 

against the longitudinal velocity fluctuation intensity, which includes 
the contributions from both turbulence and sound. A better agreement 
with the data of other investigations co\ald be obtained by plotting the 


turbulence alone. This procedure is questionable, since the sound- 
velocity fluctixations may have an appreciable effect on tiorbulent flame 
speed when coupled with tiirbiilence, and also because the magnitude of 
the sound field in the other investigations is unknown. The abscissas 
of figure 26 may be converted to tvirb\ilence intensity by the uSe of 
equation (3) and the data of table II. 

The effect of fuel-air ratio on tvirbulent flame speed of free 
flames is shown in figure 27, in which the data from the free-flame ex- 
periment reported in reference 22 are compared with those of the present 
investigation. The data from the present investigation shown in figure 
27 were taken 15.2 to 16.2 inches downstream of the 0.125-inch wire 
diameter by 0.625-inch mesh grid, while the data from reference 22 were 
taken 14 to 18 inches downstream of a 0.125-lnch wire diameter by 0.500- 
inch mesh grid. The grid sizes and axial distances are close enough 
so that the turbulence scales were roiighly equal. Because of a large 
difference in stream velocity, the two sets of data were taken under 
id.dely different velocity fluctuation intensities. Both sets of data 
show a general trend of decreasing turbulent to laminar flame-speed 
ratio in the lean and rich fuel-alr-ratio regions. The considerable 


Comparison with Data from Other Investigations 


ratio 



against the velocity fluctuation intensity 



due to 


/ 
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deviation from mean values exhibited by the data of reference 22 sub- 
stantiates the statistical deviation of turbulent flame speed noted in 
the present investigation. 


The current prevailing concept of turbulent flame propagation is 
based on references 25 and 26j wherein the proposals are made that (l) 
small-scale turbulence (of the order of the flame-front thickness) 
increases the local flame propagation rate because of increased diffu- 
sion rates^ and (2) large-scale tirrbulence increases the over-all flame 
propagation rate by wrinkling the flame with attendant increased flame- 
front area. Althou^ it is well known that a tirrbulent field contains 
a statistical distribution of "eddy sizes," the complexity of the prob- 
lem has dictated the assimiption of either the large- or small-scale case. 

The basic concept proposed by Shelkin (ref. 26) for the large-scale 
case is that the wrinkling of the flame front might be approximated by 
cones having a base proportional to the turbulence scale and a height 
proportional to the turbvilence intensity. This concept has been sup- 
ported experimentally by Hottel and coworkers in reference 1. 

Turbulent flame speeds measured in flames stabilized on bluff bodies 
in high-velocity fuel-air streams led Scurlock and coworkers to postulate 
the concept of flame -generated turbulence (refs. 6 and 7). Theoretical 
analyses to account for flame -generated turbulence were proposed by 
Karlovitz and coworkers (ref. 2) and by Sciarlock and Grover (ref. ll) , 
but neither analysis conforms to all existing experimental data. A re- 
cent analysis by Tucker (ref. 27) predicts that little, if any, turbu- 
lence would be generated by free flames with small surface deformations. 

Since the Shelkin equation is not only of historical interest but 
is also used as a basis for the Scurlock-Grover theory, the data of the 
present investigation are compared with the Shelkin equation in figure 
28, wherein the ratio S^i/Sl is plotted as a function of longitudinal 

velocity fluctuation intensity for four cumulative probabilities. The 
data include flame speeds measured with each of the three turbulence - 
generating grids and are all for the constant fuel-air weight ratio of 
0.07 (figs. 23 and 25). The Shelkin equation was correlated with the 
data at each value of cvim\ilative probability by the method of least 
squares. This correlation was accomplished by determining the constant 
B in the Shelkin equation as follows: 


Comparison of Data with Various Theories 
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The value of the constant B was estimated to be unity by Shelkin but 
varies from 0.867 for the cumulative probability of 2 percent to 0.220 
for the cumulative probability of 80 percent. The value of the con- 
stant B is shown in figure 28 for each of the cumulative probabilities. 
The data scatter about the ShelAin equation but show some agreement with 
the Shelkln theory. A curve plotted through the data obtained with the 
0.125-inch grid does not follow the shape of the Shelkln cxirve^ however. 

The theoretical expression for turbxilent flame speed presented by 
Tucker in reference 27 was obtained by a wave analysis of the effect of 
a low-intensity tiorbulent field on an inflnitesmally thin flame front 
having convolutions^ or wrinkles, of low slope. The final equation 
shows no effect of turbulence correlation coefficient and is, in the 
present notation. 


( 11 ) 

where C(t) is a function of the ratio t = Tp/Cg^. The analysis in- 
cludes the assimiption that the second-order terms (u^/s^)^ axe negll- 

J-i 

gible, that the flame tempera tvire and laminar flame speed are dependent 
only on fuel-air ratio, and that the diffusive action of turbulence is 
negligible . 

In order to obtain a comparison with the present data, the Tucker 
analysis was extrapolated well beyond its intended limit of the ratio 

(u^/Sl). The data of the present investigation compare well with the 

Tucker analysis in the low velocity fluctuation intensity region as 
shown in figure 29, but falls below the theory in the higher intensity 
range, where the theory cannot be expected to hold. A further compari- 
son of the tiirbulent flame-speed data with the Tucker analysis is made 
in figure 30, wherein the ratio S^/Sl is plotted against fuel-air 

weight ratio. These data are for each of the three tirrbulence-generating 
grids and were taken at the same mean stream velocity. The data in 
each figure fall below the theory, but, with the exception of the 
enclosed- tunnel data (fig. 30(a), fuel-air ratio, 0.045), follow the 
general trend of the Tucker analysis in the upper portion of the cumula- 
tive probability band. These data suggest that the analysis might be 
taken as the limiting value of turbulent flame speed for the turbulence 
conditions specified in the figures. 

The turbulent flame-speed data are also compared with the Scurlock- 
Grover analysis (ref. ll) in figures 29 and 30. This analysis assumes 
that the flames proceed locally at the laminar rate Sl, but that the 
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over-all flame speed &p is increased because of flame-front wrinkling 

caused by t\irb\ilence . The geometrical model is similar to Shelkin's, 
except that the wrinkle height is assumed to be positively affected by 
tiirbulent motion and negatively affected by local flame propagation. 

The analysis also accounts for flame-generated turbulence caused by 
velocity and density gradients across the flame front. 

Flame-generated turbulence has been neglected in a modified 
Scurlock-Grover analysis^ which was used in the present investigation 
and is discussed in appendix D. The measured correlation coefficients, 
as shown in figure 13, were used in the analysis, in contrast to the 
exponential forms assumed in the Scurlock-Grover report (ref. 11) and 
the comparison with free-flame data made in reference 22. The analysis 
predicted that the turbulent flame speed had not yet reached its fxill 
value at the axial stations of measurement (x = 15.2 and 16.2 in.). 

Therefore, values of Sj/Sl were calculated from the theory for flame- 

growth times corresponding to both stations and are shown in figures 29 
and 30 . • 

The data of the present investigation fall below the values of 
S|j,/Sj^ predicted by the modified Scurlock-Grover analysis but follow 

the theory fairly well in the upper portion of the probability band. 

The data follow the Sc\arlock-Grover analysis more closely for the 
largest turbulence scale (fig. 30(a)), and the agreement becomes pro- 
gressively worse as the turbulence scale becomes smaller (figs. 30(b) 
and (c)). This type of relation may be due to an increasing diffusive 
action of the turbulence as the scale becomes smaller, as discussed 
previously in this section. 

The data agree more closely with the Scurlock-Grover analysis 
than with the Tucker t-heory in the lean and rich fuel-air-ratio regions. 

This probably is due to the increasing ^/^/sL in those regions that 
tend to exceed the limits of the assumptions made in the Tucker analysis 
discussed in a preceding paragraph. 

In summary^ values from both the Tucker and the Scurlock-Grover 
analyses seem to form an upper limit to the data of the present 

investigation . 


CONCLUSIONS 

From the results of the investigation, the following conclusions 
are made : 
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1. The propagation rate of a free flame in a turbiilent field with 
intense random sound-velocity fluctuations is a statistical variable 
having a probability of occurrence band bounded by a zero propagation 
rate and possibly an upper propagation rate well above the laminar flame 
speed. 


2. The propagation rate of a free turbulent flame is substantially 
less than those previously measured for stabilized flames. 

3. The statistical spread of free-flame propagation rates is 
markedly wider in both the lean and rich fuel-air-ratlo regions. 

4. Values calculated from the Tucker analysis for free-flame propa- 
gation rates apparently form an upper limit for the statistical spread 
of free-flame propagation rates in the regions where the velocity fluc- 
tuation intensity is low. 

5. Values calc\alated from the modified Scurlock-Grover analysis 
apparently form an upper limit for the statistical spread of free-flame 
propagation rate. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland^ Ohio, April 25, 1955 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 
cross-sectional area 

smoothed time-avereige surface area of turbulent flame 
local speed of sound, ft/sec 
constant in Shelkin equation 
constant 

turbulent diffusion coefficient, sq ft/ sec 
wire diameter of turbulence-producing grid, in. 
d-c voltage, v 

rms of a-c complex wave voltaige, v 
longitudinal energy spectrum density function, sec 
longitudinal double-velocity correlation coefficient 
lateral double-velocity correlation coefficient 
hot-wire current at stream conditions, amp 
gain of decade amplifier 
longitudinal scale of turbulence, in. 
lateral scale of turbulence, in. 

Lagrangian scale of turbulence, in. 
number of flame globules 
frequency, cps 
gas constant 

Lagrangian double-velocity correlation coefficient 
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r radial distance from tunnel aucis, in. 

5 flame speed, ft/sec 

SPL sound-pressure level, db 

T temperature, °R 

t time, sec 

U mean stream velocity, ft/sec 

u longitudinal velocity fluctuation, ft/sec 

v,w lateral velocity fluctuations, ft/sec 

X fluid particle displacement due to turbulent motion, ft 

X distance downstream of turbulence-producing grid, in. 

y distance from tunnel wall, in. 

a angle, deg 

3 power band width of wave analyzer, cps 

6 flame-front thickness, ft 

5 spatial separation of velocity vectors con5)rising correlation 

coefficient, in. 

p density, Ib/cu ft 

a one standard deviation 

T ratio of flame temperature to unburned gas ten 5 )erature, Tj,/Tg^ 

(p equivalence ratio, (fuel-air ratio) /( stoichiometric fuel-air ratio) 

2 hot-wire resistance at stream conditions, ohms 

Subscripts: 

a ambient, or approach stream 

F flame 

f pertaining to correlation coefficient f 
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pertaining to correlation coefficient g 
laminar 

pertaining to a frequency interval 

no-flow condition in anemometry, or reference point in space or 
time double-velocity correlation coefficients 

pertaining to number of flames of a given radius or greater 

sound 

turbulent 

total, or time 

at a distance ? from zero position 
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APPENDIX B 


TURBULENCE AND SOUND INSTRUMENTATION AND ANALYSIS 
Theoretical Background 

Turbulent fields are generally defined by the two parameters in- 
tensity and correlation coefficient. The turbulence intensify is Refined 

as the rms value of the turbulent velocity fluctuations , and 

The double -velocity correlation coefficient may be represented by 
a second-order tensor which for a homogeneous, isotropic field of tur- 

can be reduced to two components as shown in reference 28# These 
two components are the longitudinal and lateral double-velocity correla- 
tion coefficients defined as follows: 


Longitudinal: f *= : 

u^ 


u^ 


u^ 




Lateral: g = 

V 


fo 

h- 


(Bl) 


(B2) 


where f and g are the longitudinal and the lateral correlation coef- 
ficients, respectively, u is turbulent velocity fluctuation, C is ill® 
spatial separation between the points of measuronent of the u, and the 
bar signifies a long time average. 

Longitudinal and lateral turbulence scales have been defined by the 
relations 

r^co 


Lf “ 


g 


0 


f d4 

(B3) 

g dC 

(B4) 


These scales are intended to represent a mean turbulent eddy size, and 
for homogeneous, isotropic turbulence may be related through the equa- 
tion (ref. 28) 

^ + 2 ^ ® 


(B5) 
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which, by integration, becomes 


^ 1 
2 


J 


bf 


d4 « L 


■g 


(B6) 


For the commonly observed exponential form of f ** exp/^^, equation (B6) 
■becomes ^ 


Lf = 2Lg 


(B7) 


Taylor has shown in reference 29 that the one -dimensional turbulence 
energy ^ectral density function F is the Fourier transform of the 
longitudinal correlation coefficient f: 




„ 2n:nx 
F COB — — dn 


U 


(B8) 


where F is defined by the relations 


. r^oo 

u2» u' 

Jo I 


F dn 


(B9) 


r^ 

to 

CD 

to 


and 




F dn = 1 


(BIO) 


where represents the total kinetic energy in longitudinal direction, 

and u^ represents the kinetic energy contained in the frequency interval 
n to n + dn. 


Hot-Wire Anemometry 


The constant-temperature hot-wire anemometry equipment is shown dia- 
grammatlcally in figure 31, and for the most part has been descried in 


references 16 and 30. The equation relating the rms voltage Y^t 
ured by the average square con 5 )uter to the turbulence intensity has been 
discussed in reference 16 and is as follows: 


Vi? 


4iU 


(2K) 2 (l2 - i^) 





(Bll) 
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The spectrum density function F was measured from the tape- 
recorder play-back signal by the following equation ■which is gi'ven in 
reference 16: 


The double -velocity correlation coefficient f was determined from 
a double -channel magnetic tape recording of the instantaneous voltage 
signal from the hot-wire anemometer. The double -channel tape was played 
back into a double correlation computer which electronically performed 
the instantaneous multiplication of velocity fluctuations indicated in 
equation (Bl) . The variable separation ? between velocity vectors was 
attained by the use of a movable and a fixed play-back head. The separa- 
tion 4 between the velocity vectors was equal to the separation between 
the heads multiplied by the ratio of the stream velocity to the tape 
speed. The double correlation computer was designed by D. F. Berg and 
C. C. Conger of the Le-wis laboratory. The equivalence of the correlation 
coefficient measured by this method and by the two-probe method was proven 
as a part of a research program conducted by J • C. Laurence of the Lewis 
laboratory . 


The sound instrumentation consisted of a standard microphone and 
sound-level meter from which sound-pressure level could be determined. 
The microphone signal was also impressed on the wave analyzer to obtain 
the sound kinetic energy spectrum through equation (B12) . The relation 
between sound-pressure level SPL and sound-velocity fluctuations 

has been given in reference 31 (for single harmonic waves) as 



(B12) 


where is the square of the rms voltage reading of the wave analyzer. 


Sound-Field Instrumentation 



(B13) 


32 


NACA TN 3456 




APPENDIX C 


PHOTOMULTIPLIER INSTRUMENTATION 


Relation Between Flame Light Intensity and Flame Speed 


Assumptions. - The conditions assumed in relating flame light in- 
ten si tyliljSrcated hy the voltage output of a photomultiplier tube to 
the surface area of a flame axe treated in reference 32 . They may be 
summaxized as follows: 


(1) It is postulated that, if a point source of light is increased 
to a finite size but the distance from the light source to the detector 
is constant and large in comparison with the maximum dimension of the 
flame, the light flux inpinging on the detector will obey the inverse 
square law. Under this condition, a change in the surface area of the 
flame will produce an approximately linear change in the flux intensity 
registered by the detector . 

(2) The flame must consist of a homogeneous cloud of nonabsorbing 
emitters. 

If the propagation of a free-flame globule moving in a gas stream at 
stream velocity is observed by a detector through a long slit placed 
axially with respect to globule travel, some additional requirements are: 

(3) The detector has an essentially flat light response over the 
entire field of view . 


(4) The slit is narrow enough in relation to flame -globule diameter 
to approximate a diameter of the projected spheroid. 

(5) The flame-globule center must not suffer any appreciable radial 
di^lacement from the tunnel centerline. 

(6) The flame globule must not have a preferential growth direction. 


(7) The flame temperature must be constant for all the flame 
globule s . 

Conditions (l) and (4) were met by suitable placement and dimension- 
ing of the apparatus. For condition (3), it was determined that the 
photomultiplier unit had a maximum light response error of 5 percent over 
the total light acceptance angle. Photographic evidence obtaine n e 
preliminary phase of the investigation Indicated satisfactoiy fulfillment 
of requirements (5) and (6) . Assumption (7) has been previously treated 
in the section of the DISCUSSION entitled Measurements of Present 
Inve sti gat ion • 
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Description of apparatus . - A 931A type photomultiplier having an 
adjustable-plate voltage power supply was mounted in a light-tight box 
behind a 0.25-inch-dlameter aperture. This detector unit was positioned 
in an enclosed tunnel 4 feet from the free-jet stub tunnel, as shown in 
figure 5 • The detector aperture was centrally located with respect to 
the slit in the machined aluminum plate, which replaced an upper portion 
of one wall of the stub tunnel. The voltage output of the photomultiplier 
was an 5 )lified and fed to an oscilloscope. Oscillograms were teiken by 
high-speed, continuous motion picture technique. 

Oscillograms . - As shown in figure 20, the oscillograph trace of the 
ignition and growth of a single flame consists of two portions, a spark 
light amplitude decay curve and a flame light amplitude growth curve 
whose slope at any point is equal to the change of diameter of the flame 
globule with time, governed by the considerations listed previously. 

This slope 2dr/dt is proportional to turbulent flame speed Sij as 
given by equation (2) . 


Relation Between Spark Light Intensity and Spark Current 

< 

o , , 

Standring and Looms (ref. 33), using short-duration sparks produced 

by a lightly damped oscillatory current, have shown that the light emitted 
by a given spark gap depends on the current in the gap and its time 
vaxlation . 

In order to determine the effect of a long-duration, direct-current 
pulsed spark on the current - spark light relation, an auxiliary investi- 
gation was undertaken in collaboration with Clyde C. Swett, Jr. The 
spark production and recording apparatus is fully described in reference 
34. The photomultiplier detector apparatus was as described previously 
except that the detector unit was mounted at one end of a rectangular 
tube 2 feet in length. The field of view of the photomultiplier tube 
was limited by a 0.0625- by 2-inch aperture in the other end of the 
tube . A Polaroid Land camera was used to make oscillograms . 

A series of oscillograms was taken of ^arks ignited in still air 
and in an air stream having a velocity of 50 feet per second. Spark 
duration was about 600 microseconds, spark energy was of the order of 
17 mllli joules, and peak spark current approximated 0.23 amperes. A 
comparison of spark light amplitude -time curves with the current-time 
curves obtained in the apparatus designed by Swett showed a marked simi- 
larity of shape between the sets . It was concluded from this information 
that, at least for a first-order approximation, spark light intensity was 
roughly proportional to spark current . 
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APPENDIX D 


CALCULATION OF FLAME SPEED FROM SCURLOCK-GROVER ANALYSIS 


The basic equation in the Scurlock-Grover theory (ref. Il) is ob- 
tained by assinoing that the increase of flame speed is due to the wrin- 
kling of the flame front in the form of right circular cones, so that 


Sip 




(Dl) 


where ^ /Lg is assumed proportional to the ratio of cone height to 
base diameter and C]^ is a constant taken equal to 4. The mean square 
particle displacement is assumed to depend on three effects: 

(1) Tixrbulent mixing, which increases the flame-front area 

(2) T.aTTi1rmr flame propagation, which decreases the flame-front area 

(3) Flame -generated turbulence, which increases the turbulent mixing 
The first effect on X^ is given by 



(D2) 


where ^ is the Lagrangian correlation coefficient and g^ is the 
Eulerian correlation coefficient. The second effect is given by 


where C2 
first two 


dt 


= Cg Sl' 


IS a constant t 
effects is then 


_ 


dt 


( / 9 \ 

1/2 ' 

' ' (3) 

- 1 ^ 

\}- 

> 

equal to l/4. 

The coml: 

^ d x| ( 


2'\j4 

dt 
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Qince "the Lagranglan corrcXation coefficienb canno't ba measurad. d-irecbly^ 
Scurlock and Grover assumed it to be of the form given in reference 15: 


vhere ^ is the Lagrangian scale of turbulence and was assumed to be 
equal to l/2 Lg. In the present analysis, 31 was assmned^o be equal to 

f, through the independent variable transformation K =Vn2 t. Isotropy 
■was assumed^ so ‘bha'b g was olD'tained from f ‘through 'the ralafion (rof • 

28) 


The congelation coefficient was then obtained from the independent 

variable transformation 4 = S^t. The tiirbulence scale Lg was calcu- 
lated by matching the curve exp(4/Lg) to the g plot previously calcu- 
lated from f by equation (D6) . 

All integrations and differentiations were done graphically, and the 
calculation procedure followed the following order: 

(l) Calculate g from equation (D6). 



(D5) 


(2) Transform 

(3) Transform g(4) bo . 



(4) Calculate 



(5) Calculate 



1/2 


(6) Calculate 
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(7) Calculate 



(8) Calculate 



(9) Calculate 


1 



As mentioned in the DISCUSSION section^ flame-generated turbulence 
has been assumed negligible, and the final form of the Scurlock-Grover 
theoretical equation was as given above. 
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TABLE I. - SUMMARY OF LONGITUDINAL TURBULENCE SCALES 

[Estimated from longitudinal energy spectra in enclosed- 
tunnel installation. Mean stream velocity, 50 ft/sec; 
static pressure, 26 in. Hg abs; static temperature, 540° 
R. Comparative values calculated through the isotropic 
turbulence relation L^ = 2Lg from lateral scale of t\ir- 

bulence data reported in ref. 15.] 


Distance downstream 
of grid, X, 

in. 

Longitudinal scale of turbulence^ 
in. 

Present in- 
vestigation 

Ref. 15 

6.3 

0.16 

— 

9.3 

.21 

— 

12.3 

.25 

0.17 

15.3 

.25 

.21 
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TABLE II. - COMPARISON OF VELOCITY FLUCTUATIONS IN TURBULENT STREAM 
WITH AND WITHOUT INTENSE SOUND FIELD 


[Velocity fluctuations measured along duct axis. Measured value of 
sound -pres sure level, 136 db; stream velocity, 68.5 ft/sec. j 


Grid 

size^ 

in. 

Distance down- 
stream of 
turbulence- 
producing 
grid^ X, 
in. 

stream velocii: 
ft/ 

iy fluctuations, 
^sec 

Calculated values of 
sound field 

Measured values of 
sound field 

Diam. 

Mesh. 

With sound 
field, 

vS 

Without sound 
field, 

Vi 

Vi = V u2 . 

ft/ sec 

Sound- 

pressure 

level, 

db 

Vi, 

ft/sec 

Sound- 

pressure 

level, 

db 

0.125 

0.625 

3 

A 

1.81 

1.40 

1.15 

137.5 

0.96- 

136 



4: 

1.75 

1.32 

1.15 

137.5 

.96 

136 

0.063 

0.313 

3 

A 

1.32 

0.926 

0.94 

135.5 

0.96 

136 



4 

1.41 

.892 

1.10 

137.0 

.96 

136 

0.0313 

0.156 

3 

A 

1.12 

0.688 

0.89 

135.0 

0.96 

136 



4 

1 . 22 

.619 

1.05 

136.5 

.96 

136 
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(a) Enclosed-tunnel installation. 

Figure 1. - Air-ducting, wind-tunnel, and mean air-flow instrumentation. 
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(b) Free -jet installation. 

Figure 1. - Concluded. Air-ducting, wind-tunnel, and mean air-flow instrumentation. 
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Figirre 2. - Fuel system. 
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Figure 3. - Photographic instrumentation for measurement of free-flame growth in 
enclosed-tunnel installation. 
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(a) Over-all view of installation- 

Figure 4. - Ionization- gap instrumentation for measurement of free- 
flame growth in free- jet installation. 
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(b) Instrumentation 

Figure 4. - Continued. Ionization- gap instrumentation for measurement of 
free -flame growth in free -jet installation. 
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(c) Ionization-gap probe. 



(d) Construction details of ionization-gap probe. 

Figure 4. - Concluded. Ionization-gap instrumentation for measurement of 
free-flarae growth in free- jet installation. 
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Figure 5. - Photonrultiplier instrumentation Tor simultaneous measurement of free- 
flame growth and spark light intensity. 
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Tunnel traverse, in. 


(a) Enclosed-tunnel installation. Nfeasured 10.75 
inches downstream of grid of 0.125-inch-diameter 
wire by 0.625 -inch mesh. 

Figure 6. - Typical mean stream velocity profiles. 
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(t>) Free-Jet installation. Reference velocity equals velocity at y = 0.5 
Inch and x = 1.0 inch downstream of spark electrodes. 


Figure 6. - Concluded. Typical mean stream velocity profiles 
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Figure 7. - Typical fuel -air -ratio profiles in free- jet installation. 
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Figure 8. - Longitudinal turbulence intensity measiared in enclosed- tunnel Installation. 
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Longitudinal energy spectrum density function, F, sec 
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Figure 9. - Typical longitudinal energy spectrum of turbxilence in enclosed- tunnel 
installation. Mean stream velocity, 50 feet per second. 
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(a) Grid: 0.125- inch wire diameter, 0.625- inch mesh. 

Figure 10. - Longitudinal velocity fluctuation intensity in free-Jet installation. 
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(b) Grid: 0.063- inch wire diameter, 0.313- inch mesh. 


Figure 10. - Continued. Longitudinal velocity fluctuation intensity in 
free- Jet installation. 
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(c) Grid: 0.0313- inch wire diameter, 0.156- inch mesh. 


Figure 10. - Concluded. Longitudinal velocity fluctuation intensity in 
free- jet installation. 
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Figure 11. - Longitudinal kinetic- energy spectrum of free- jet sound field. 
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Figure 12. - Typical longitudinal velocity fluctuation energy spectriim measured in free- Jet 
installation. Grid; 0.125- inch wire diameter, 0.625- inch mesh; distance downstream of 
spark electrodes, 4 inches; mean stream velocity, 68.5 feet per second; static pressure, 
29 Inches of mercury absolute; static temperature, 545° R. 
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Figure 13. - Longitudinal double-velocity correlation coefficients measured in free- Jet installation. 
Mean stream velocity, 68.5 feet per second. 
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(b) Grid; 0.065-lnoh wire diameter, 0.313-lnch mesh. 


Fi^re 13. - Continued. Longitudinal double-velocity correlation coefficients measued 
Installation, Mean stream velocity, 68.5 feet per second. 
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(c) Grid: 0. 0313-Inch wire diameter, 0.153-lnch mesh. 

Figure 13. - Concluded. Longitudinal double-velocity correlation coefficients measured In free-Jet 
installation. Mean stream velocity, 68.5 feet per second. 
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(a) Typical photographic record. Film speed. approximately 3600 
frames per second. 

Figure 14. - Typical turbulent flame- speed data from enclosed- tunnel 
installation. 
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(c) Cumulative probability of turbulent flame speed. Mean stream velocity, 60 
feet per second; propane-air ratio, approximately 0.045 by weight; static 
pressure, 26 Inches of mercury absolute; temperature, 540° R. 


Figure 14. - Concluded. Typical turbulent flame-speed data from enclosed- 
tunnel installation. 
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(a) Fuel-air ratio^ 0.056; mean stream velocity, 68.3 feet 
per second. 


— Figure 16. - Typical cumulative probability of flame-globule size 
measured vith ionization-gap instrumentation in free- Jet in- 
stallation. Grid: 0.125- inch wire diameter, 0.625- inch mesh; 

static pressure, 29.3 inches of mercury absolute; static tem- 
perature, 545° R. 
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(b) Fuel-air ratio^ 0.070; mean stream velocity^ 69.6 feet 
per second. 

Figure 16. - Concluded. Typical cumulative probability of flame- 
globule size measured with ionization-gap instrumentation in 
free- jet installation. Grid: 0.125- inch wire diameter, 0.625- 

inch mesh; static pressure, 29.3 inches of mercury absolute; 
static temperature, 545® R. 
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Figure 17. - Typical flame-globule diameter histories measured In enclosed- 
tunnel installation. 
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Figure 18. - Typical oscillograph of ionization -gap and hot-wire signals in free- jet stream at time of flame-glohule passage. 
Grid: 0.125-inch wire diameter, 0.625-lnch mesh; stream velocity, 70 feet per second; static pressure, 29 Inches of 

mercury absolute; temperature, 80° F. 
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Figure 19. - Effect of spark source energy on flame -globule 
growth. 
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Figure 20. - Typical oscillographs of photomultiplier prohe trace depicting spark occurrence and growth of flame 
globule in free -jet stream. 
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Figure 21. - Correlation between spark light Intensity and turbulent flame speed for 
151 consecutive flame globules at identical mean stream and fuel conditions. 
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Figure 22. - Statistical variation of flame 
speed as determined from photomultiplier 
oscillograms . 
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Figure 23. - Turbulent flame speed as a function of longitudinal velocity fluctuation intensity with cumulative 
probability of occurrence as a parameter. Propane-air weight ratio, 0.070; stream static temperature, 545^ P; 
static pressure, 29.3 inches of mercury absolute; grid: 0.125- inch wire diameter, 0.625-inch mesh; free- Jet 

installation. 
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Figure 24. - Turbulent flame speed as a function of turbulence Intensity. Propane-air weight 
ratio, approximately 0.045; stream temperature, 540*^ R; static pressure, 26 Inches of mer- 
cury absolute; enclosed-tunnel Installation; grid; 0.125-lnch wire diameter, 0.625-lnch 
mesh . 
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Figure 25, - Turbulent flame speed as function of fuel-air ratio with cumulative 
probability of occurrence as parameter. Mean stream velocity, 70 feet per 
second; stream static pressure, 1 atmosphere; static temperature, 85° F. 
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(b) Grid; 0.063-inch diameter, 0.313-inch mesh. 

Figure 25. - Continued. Turbulent flame speed as function of fuel-air 
ratio vrith cumulative probability of occxirrence as parameter. Mean 
stream velocity, 70 feet per second; stream static pressure, 1 at- 
mosphere; static temperature, 85° F. 
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Fuel-air weight ratio 

(c) Grid: 0.0313-inch diameter, 0.156-inch mesh. 

Figure 25. - Concluded. Turbulent flame speed as function of fuel-air 
ratio with cumulative probability of occurrence as parameter. Mean 
stream velocity, 70 feet per second; stream static pressure, 1 at- 
mosphere; static temperature, 85^ F. 
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(a) Free-jet data. Propane-air ratio, 0.070. 

Figure 26. - Cotiiparison of turbulent flame-speed data with those of other 
investigations . 
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(b) Enclosed- tionnel data. Propane- air ratio, approximately 0.045. 

Figure 26. - Concluded. Con5)arison of turbulent flame-speed data with those 
of other investigations. 
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Figure 27. - Comparison of turbulent flame-speed data with those of reference 22 
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(a) Cumulative probability of occurrence, 2 percent. 
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(b) Cumulative probability of occurrence, 20 percent. 

Figure 28. - Comparison of turbulent flame -speed data with 
theory of reference 26. Fuel-air weight ratio, 0.07. 
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(d) Cumulative probability of occurrence, 80 percent. 

Figure 28. - Concluded. Comparison of turbulent flame-speed data 
with theory of reference 26. Fuel-air weight ratio, 0.07. 
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Figure 29. - Comparison of turbulent flame-speed data with theories of references 11 and 27. Propane-air weight 
ratio, 0.07; stream static pressure, 29.3 inches of mercury absolute; temperature, 545° R; grid; 0.125- inch 
wire diameter, 0.625-inch mesh. 
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(a) Grid: 0.125-inch wire diameter, 0.625-inch mesh. 

Figure 30. - Comparison of turbulent flame-speed data with theories from references 11 and 27. 
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(b) Grid: 0.063-inch wire diameter, 0.313-inch mesh. 


Figure 30. - Continued, 
references 11 and 27. 
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(c) Grid: 0.0313- inch wire diameter, 0.156- inch mesh. 

Figure 30. - Concluded. Comparison of turbulent flame speed data with theories from 
references 11 and 27. 
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